o 



43 



o 



>< 



HYDRODYNAMIC LIMIT FOR A TYPE OF EXCLUSION 
PROCESSES WITH SLOW BONDS IN DIMENSION > 2 

TERTULIANO FRANCO, ADRIANA NEUMANN, AND GLAUCO VALLE 



Abstract. Let A be a connected closed region with smooth boundary con- 
tained in the d-dimensional continuous torus T d . In the discrete torus N~ 1 TT^ r , 



00 _ we consider a nearest neighbor symmetric exclusion process where occupancies 

of neighboring sites are exchanged at rates depending on A in the following 
way: if both sites are in A or A , the exchange rate is one; If one site is in A 
and the other one is in A^ and the direction of the bond connecting the sites 
C^ is Bj, then the exchange rate is defined as A -1 times the absolute value of the 

CLj ■ inner product between ej and the normal exterior vector to dA. We show that 

this exclusion type process has a non-trivial hydrodynamical behavior under 
diffusive scaling and, in the continuum limit, particles are not blocked or re- 
J^ , fleeted by dA. Thus the model represents a system of particles under hard 

core interaction in the presence of a permeable membrane which slows down 
the passage of particles between two complementar regions. 



1. Introduction 

> 

^^ ' The exclusion process is a continuous time interacting particle system where 

f*"«« , particles move as independent random walks on a graph except for the exclusion 

rule that prevents two particles from occupying the same site, or vertex. In the 
symmetric case, the process evolves as follows: to each bond we associate a waiting 
exponential time, which are independent of the waiting time for any other bond; 

^^ , at the waiting time the occupancies of the sites connected by the bond are ex- 

changed; the parameter of the exchange times, or exchange rate, depends only on 
the bond. The especification of the exchange rates determines the environments 
for the exclusion process. In our case, as the underlying graph, we consider the 
discrete torus with N d points and nearest neighbor bonds. The variable N is the 
5^ ■ scaling parameter. 

This paper studies the hydrodynamical behavior of symmetric exclusion pro- 
cesses in non-homogeneous environments with the presence of slow bonds. Here 
we mean non-constant environments where a usual bond has exchange rate one 
and a slow bond has exchange rate lower than one. With respect to the scaling 
parameter, we assume that a slow bond has exchange rate of order N^ 1 . When the 
environment is constant, the exclusion process has a well-known hydrodynamical 
behavior under diffusive scaling, but, in the presence of slow bonds, particles will 
not move fast enough to garantee that we still have an hydrodynamic behavior in 
diffusive scaling. Hydrodynamics in diffusive scaling have been obtained in several 
cases, even when the environment is random and consists only of slow bonds. 

For one dimensional processes, in [5], the exchange rate over a bond [Jr, ^=] is 
given by [N(W(x + 1/N) — W (x / N))]~ l : where W is an a-stablc subordinator of 
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a Levy Process. They obtain a quenched hydrodynamic limit. In papers previous 
to [5], for example |3] and [9], the randomness or non-homogenity did not survive 
in the continuum limit. An also one-dimensional result, following [5], was obtained 
in [6], for more general, but non-random, increasing functions W. The techniques 
used in those papers were strongly based on theorems about convergence of one 
dimensional continuous time stochastic processes. In fact, even the d-dimensional 
case treated in [10J has considered a class of non-homogeneous environments that 
could be decomposed, in a proper sense, in d one-dimensional cases. Recently, 
different approaches have been searched to deal with d-dimensional environments, 
see [4] and [7]. 

We now describe the exclusion processes we are concerned. Let {ej : j = 1, ..., d} 
be the canonical base of R d and A C T d be a simple connected region with smooth 
boundary dA. If the bond [-^, —jp-] e A _1 T^ has vertices in each of the regions 
A and A , its exchange rate is defined as A^ -1 times the absolute value of the inner 
product between ej and the normal exterior vector to dA. For others edges, the 
exchange rate is defined as one. This means that the slow bonds are among those 
crossing the boundary of A. We call this process the exclusion process with slow 
bonds over dA. 

Wc can interpret <9A as a permeable membrane, which slows down the passage 
of particles between the regions A and A . For this type of exclusion process, the 
membrane does not completely prevent the passage of particles, and still survives 
in the continuum limit, appearing explicitely in the hydrodynamic equation. The 
exchange rate of particles for a bond crossing dA is smaller if the bond is close 
to a tangent line of dA. Note that this assumption has physical meaning, take 
for example cases of reflections in several physical models: partial reflection of 
light crossing a media with diferent refraction indexes, mechanical systems where 
particles try to cross some interface, etc. However the direction of the speed of 
particles is not changed as usually occur in physical reflection. Our definition of 
the exchange rates also allows a strong convergence result for the empirical measures 
associated to the exclusion process making simpler the proof of the hydrodynamic 
limit. 

The hydrodynamical equation of the exclusion process with slow bonds over dA is 
a parabolic partial differential equation dtp = Cap, where the operator C\ is a sort 
of d-dimensional Krein-Feller operator. Without the presence of slow bonds, the 
operator £a would be replaced by the laplacian operator acting on C 2 functions and 
the hydrodynamical equation is therefore the heat equation. Here, the existence 
of the membrane modifies the domain, and thus the operator itself. In fact, we 
observe that the proper domain for £j\ contains functions that are discontinuous 
over dA. Geometrically, C\ glues the discontinuity of a function around dA and 
then behaves like the laplacian. 

One possible approach to prove the hydrodynamic limit for the exclusion process 
with slow bonds over dA is through Gamma convergence. In [7] , this approach and 
the conditions for it to hold are discussed, see also [3]. There, the coersiveness 
condition would require some kind of Rellich-Kondrachov's Theorem (namely, the 
compact embedding in L 2 of some sort of Sobolev space supporting an extension 
of C\, see [2]). In the method presented here, we go in this direction, but instead 
of reach the hypotheses in [7] , we have use similar the analytical tools to obtain a 
short and simple proof of uniqueness of the hydrodynamic equation. Wc also show 
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that the extension of C\ satisfies the Hille-Yoshida Theorem. On the other hand, 
the convergence from discrete to continuous that we present here is made in very 
direct way, and it was inspired by the convergence of the discrete laplacian to the 
continuous laplacian. 

The paper is presented as follows: In Section [3J we define the model and state 
all results proved to be proved in the paper; Section [3] is devoted to prove all the 
results concerning to the continuous operator Ca', In Section 21 the hydrodynamic 
limit is proved. 

2. Notation and Results 

Let T d be the d-dimensional torus, which is [0, l) d with periodic boundary con- 
ditions, and T^ be the discrete torus with N d points, which is {0, ..., N — l} d with 
periodic boundary conditions. We denote by r\ = (rj(x)) x£ jd a typical configuration 

in the state space J7jv = {0, 1} T ™, for which, r](x) = means that site x is vacant, 
and rj(x) = 1 that site x is occupied. If a bond of iV -1 T^ has vertices fj and 44, 
it will be denoted by [■&, 44]. 

Recall that {ej : j = 1, ..., d} is the canonical base of K d . The symmetric nearest 
neighbor exclusion process with exchange rates ^ y > 0, x, y G T 1 ^, \x — y\ = 1, 
is a Markov Process with configuration space VLn, whose generator Ljv acts on 
functions / : 51 n ->Ras 

d 

(L N f)(v) = E E C+* [f(v x>x+ej ) - m] , (2.i) 

x£T% 3 = 1 

where r] x,x+ej is the configuration obtained from r\ by exchanging the variables r\{x) 
and r\(x + ej): 



r](x + ej), if y = x . 
(v x ' x+ej )(y) = < V(x), Hy = x 



e 



j ■■ 



rj(y), otherwise. 

Let v£ , a £ (0, 1), be the Bernoulli product measure f2jv ; i-e., the product measure 
whose marginals have Bernoulli distribution with parameter a. Then {i/^ ■ < 
a < 1} is a family of invariant, in fact reversible, measures for any symmetric 
exclusion process. 

Now, fix a simple connected region A C T d with smooth boundary dA. Denote 
by C( u ) the normal unitary exterior vector to the smooth surface dA in the point 
ii£9A. Iff £ A and 2+£i e A c , or § G A c and -^- G A, we define ( x ,j 
as a vector C( u ) evaluated in an arbitrary but fixed point u G dA n [x, x + ej]. 
The exclusion process with slow bonds over dA is a symmetric nearest neighbor 
exclusion process with exchange rates £xx+e- ~ £x+ex gi ven by 

^_M if » e A and -^ G A c , or f G A c and -p- G A, 



iV 
1 , otherwise 



(2.2) 



for j = 1, . . . , d, and for every x G Tjy. In this case, the exchange rate of a bond 
crossing the boundary dA is also of order A^ 1 , but it depends on the angle of 
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Figure 1. The darker region corresponds to A. The bolded bonds 



have exchanges rates 
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- , any other bond has exchange rate 1 . 



incidence: the crossing of <9A by a particle gets harder to happen as the direction 
of entrance gets closer to the tangent plane to the surface dA. 

From now on, the rates in the definition of Lm will always be given by 12.21 
Denote by {rj^ : t > 0} a Markov process with state space fijv and generator L^ 
speeded up by A^ 2 . Let £)(R + ,fijv) be the Skorohod space of cadlag trajectories 
taking values in $7 at. For a measure fi on fi;v, denote by P^ the probability measure 
on D(R+, Oat) induced by the initial state \i and the Markov process {r)f : t > 0}. 
The expectation with respect to P^ is going to be denoted by Ej*. 

A sequence of probability measures {fi^ : N > 1} is said to be associated to a 
profile 7 : T d —> [0, 1] if /ijv is a probability measure on Q N , for every N, and 



lim 

N—KX 



(J-N 






H(u)~/(u)du 



>sy = o 



(2.3) 



for every S > 0, and every continuous functions H : T d — » R. 

The exclusion process with slow bonds over 9A has a related random walk on 
A^ _1 T A > that describes the evolution of the system with a single particle. Thus 
particles in the exclusion process evolve independently as such random walk except 
for the hard core interaction. To simplify notation later, we introduce here the 
generator of this random walk, which is given by 



^NH)(§)=J2{C, X+ e 3 [tf(^i) "#(#)] +€^_ ej [H(^l)-H(§)]}, 



N- X T% 



for every H 

for functions H defined on T d and on A" 



and every x £ Tjj . We will not differentiate the notation 



■N 
-ljd 



N- 



2.1. The Operator C\. Here we define the operator £\ and state its main proper- 
ties. Its domain is defined as a set of functions two times continuously differentiable 
inside and outside A, but possibly discontinuous on dA. Besides, in the boundary 
d A, wc should require particular conditions of those functions in order to have 
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good properties of £a that allows us to conclude the uniqueness of solutions of the 
hydrodynamic equation and obtain a strong convergence result for the empirical 
measures in the proof of the hydrodynamic limit. The necessity of these conditions 
are going to be made clear later in the text. 

Definition 2.1. Recall that £ denotes the normal exterior vector to the surface 
dA. The domain £> A C L 2 (T d ) will be the set of functions H G L 2 (T d ), such that 
H(u) = h(u) + A1a(u), where: 
(i) A G K; 
(ii) h G C 2 (T d ); 
(iii) Vh\ dA {u) = -\({u). 
Now, we define the operator C\ : 3D a —> L 2 (T d ) by 

C A H = A/i . 

Geometrically, the operator C\ removes the discontinuity around the surface d A 
and then acts like the laplacian operator. 

Remark 2.1. It is not entirely obvious why there exist functions h G C 2 (T d ) such 
that S/h\d\(u) = —\£(u), for A ^ 0. For an example of such a function, consider 
firstly g : T d — ► R defined by 

Xdist(u,dA) , if u G A c , 



9 ^ \ -Xdist(u,dA), if ueA. 

One can checks that this function is smooth in an open neighborhood V of dA, and 
satisfies the condition \7g\g\(u) = — A£(u). However, g is not differ entiable in 
the space T d . To solve this problem, it is enough to multiply g by £^ $,, where 
{$i} is a partition of unity such that the support of any $i is contained in V and 
J2i^i( u ) ~ 1 f or all u ^ U C V , U an open set containing dA. Finally, the 
function 






h{u) = g(u) > $i(u) 






satisfies the required conditions. 



For the next result we need to introduce some notation. We denote by I the 
identity operator in L 2 (T d ) and by ((•, •}) and || • || its usual inner product and norm: 

«/,.9» = / f(u)g(u)du and ||/|| = VT7J)) , f, 9^L 2 (T d ). 

Theorem 2.2. There exists a Hilbert Space {T-L\, ((•, -))i,a) compactly embedded in 
L 2 (T d ) such that S)a c H\ and £a can be extended to C\ : H\ — > L 2 (T d ) in such 
way that the extension enjoys the following properties: 

(a) The domain T-i 1 ^ is dense in L 2 (T d ); 

(b) The operator C\ is self-adjoint and non-positive: ((H, —C^H)) > 0, for all 
H in H\; 

(c) The operator I — C\ : H\ — > L 2 (T ) is bijective and 2) a is a core for it; 

(d) The operator C a is dissipative, i.e., 

\\liH-£ A H\\>n\\H\\, 
for all H G H\ and jjl > 0; 
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(c) The eigenvalues of —£\ form a countable set = po < pb\ < • • • with 

lim n _s.oo p n = oo, and all these eigenvalues have finite multiplicity; 
(f) There exists a complete orthonormal base of L 2 (T d ) composed of eigenvec- 
tors of —Ch- 
in view of (a), (c) and (d), by the Hille- Yoshida theorem, C\ is the generator of a 
strongly continuous contraction semigroup in L 2 (T d ). 

The notation TL\ will be defined in Section [31 and it has been chosen in analogy 
to the notation used for Sobolev spaces. 

2.2. The hydrodynamic equation. Consider a bounded Borel measurable profile 
Po : T d — > K. A bounded function p : M + x T d — > M is said to be a weak solution of 
the parabolic differential equation 

d t p = C A p 



p(o,0 = poO), (2 ' 4) 

if for all functions H in H\ and all t > 0, p satisfies the integral equation 

((p t ,H)) - ((p ,H))- f {{ Ps ,C K H))ds = 0, (2.5) 

where pt is the notation for p(t, •). We prove in Subsection 14.31 the uniqueness 
of weak solutions of 14.271 Existence follows from the convergence result for the 
empirical measures associated to the diffusively scaled exclusion processes with 
slow bonds over A, this is discussed in Section [4] 

Theorem 2.3. Fix a Borel measurable initial profile 7 : T d — >• [0, 1] and consider 
a sequence of probability measures /ijv on Qn associated to 7. Then, for any t > 0, 



JknJ^II^ J^ H{x/N) m {x) - J H(u)p(t,u)du\> s] 



xGTi 



= 0. 



for every 5 > and every function H £ C(T d ), where p is the unique weak solution 
of the differential equation (|2.4[) with po = 7- 

3. The operator £a 

We begin by studying properties of £a defined in the domain T>\ and we consider 
the extension afterwards. 

Lemma 3.1. The domain 35a is dense in L 2 (T d ). 

Proof. It is enough to prove that there exists a subset of 33 a which is dense in 
L 2 (T d ). All smooth functions with support contained in T d \i9A belong to 33a, 
which is clearly a dense subset of L 2 (T d ), since <9A is a smooth zero Lebcsguc 
measure surface that divides T d \9A in two disjoint open regions. □ 

From now on, we use £4 to denote the d-dimensional Lebesgue measure on T d . 

Lemma 3.2. The operator —C\ : 33a — > L 2 (T d ) is symmetric and non-negative. 
Futhermore, it satisfies a Poincare inequality, which means that there exists a finite 
constant C > such that 

\\H\\ 2 <C{{-C A H,H)) + n H(x)d x y (3.6) 
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for all functions H £ 'S A . 

Proof. Let H,G E 1) A . Write H = h + Xh 1a and G = g + X g 1a, as in definition 
12.11 By the first Green identity and condition (iii) in definition 12. 1[ we have that 

X h [ Agdu = X h [ {Vg-C)dS = -X h X g Vol rf _i(dA) (3.7) 

J A JdA 

= X g I (V/i • C) dS = X g I Ah du , 

JdA J A 

where dS is a infinitesimal element of volume of <9A and Vo\d-i(dk) is its (d — 1)- 
dimensional volume. Thus, 

((H, -C A G)) = ((h + X h l A , -Ag)) = - f hAgdu-X h f Ag du 

JT d J A 

= - [ gAhdu-X g J Ahdu = ({-£ A H,G}). 

JT d J A 

For the non-negativeness, using l3~7l above. 

((H,-£ A H)) = - h Ahdu - X h f Ahdu 

JT d J A 

|V/i| 2 du + A^Vol d _i(9A) > 0. 



It, i 
It remains to prove the Poincare inequality. Write 

ll#H 2 -( / H(x)dx) = [ \h(u)- [ H(v)dv 
\Jjd I Jfd L J T d 

which can be rewritten as 

(h{u) - J ^ h(v) dr^j + X h (l A (w) - 4(A))] 

Now apply the inequality (a + b) 2 < 2 (a 2 + b 2 ) to the previous expression to obtain 
that it is bounded by 

2 / (h(u)- f h(v)dv) 2 du + 2X 2 U d {K)-{(i d (K)) 2 ). 

Jjd \ Jjd / v / 

By the usual Poincare inequality, see [2] , the last expression is less than or equal to 



j, i 



du . 



du . 



2d / \Vh(u)\ 2 du + 2X 2 h (t d (A)-(t d (A)y 



lid 

Choosing a constant C 2 > such that £ d {A) - (4(A)) 2 < C 2 VoLi-i(<9A), the 
previous expression is bounded above by 



2max{C 1 ,C 2 }«-£ A #,tf», 
which finishes the proof with C = 2 max{Ci, C2}. □ 

Denote by ((■, -))i.a the inner product on 2Da defined by 

((F,G))i, A = <(F,G» + ((F,-C A G)). 

Let 7i A be the set of all functions F in L 2 (T d ) for which there exists a sequence 
{in : n > 1} in Da such that F n converges to F in L 2 (T d ) and F n is Cauchy for 
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the inner product ((•, -))i,a- Such sequence {F n } is called admissible for F. For F, 
G in H\, define 

{{F,G)) 1 . A = lim((F„,G„»i, A) (3.8) 

n— >oo 

where {F n }, {G n } are admissible sequences for F, G, respectively. By [TTJ Proposi- 
tion 5.3.3], the limit exists and does not depend on the admissible sequence chosen. 
Moreover. H\ endowed with the scalar product ((•, -))i,a just defined is a real Hilbcrt 
space. From now on, we consider 7i\ with the norm induced by ((■, -}}i,a, unless we 
mention that we are going to use the L -norm. 

Lemma 3.3. The embedding H\ C L 2 (T d ) is compact. 

Proof. Let {H n } a bounded sequence in %\. Fix {F n } as a sequence in 15a such 
that \\F n — H n \\ — >• and {F n } is also bounded in T~L\. Thus, to get a convergent 
subsequence of {H n }, it is sufficient to find a convergent subsequence of {F n } in 
L 2 (T d ). Write F n = /„ + A„1 A , with /„ e C 2 (T d ). Then, 

((■Fn, F„))i, A = ((/„ + A„1 A , /„ + A„1a)) + ((/„ + A„1a, -A/ n » . 
Expanding the right hand side and using [3~7"1 we get that 

((F n ,F n )) 1A = \\f n \\ 2 + A 2 4(A) + 2A„ [ /„(«) du + || V/„|| 2 + A 2 Vol^^aA) , 

JA 

which is greater or equal to 

||/„|| 2 + A 2 ^(A) - A 2 - 4(A) / fl(u) du + \\V.f n f + A 2 VoLj-i^A) 

JA 

= (^(A)-l + Vol d -i(aA))A 2 +(l-£ d (A))y f 2 (u)du + Jj 2 (u)du+\\Vf n \\ 2 

> (Vold_i(5A) - £ d (A c )) A 2 + (1 - £ d (A)) ||/ n || 2 + ||V/„|| 2 . 

If we put /„ = /„ + 1, and write F n = /„ — A„l A c , an analogous computation shows 
that ((f„,F„))i : A is greater or equal than 

Vol d _i(dA) - 4(A)) A 2 + (1 - 4(A C )) ||/„|| 2 + ||V/„|| 2 . 

By the classical isoperimetric inequality on the Torus (see [1] Lemma 4.6] for the 
statement and a direct proof), we have that 

max{ VoLi_i(c>A) - 4(A C ) , Vohi_i(dA) - 4(A) } > . 

Since {{{F n , F n ))i t \} is a bounded sequence, we conclude that {A„} is bounded, 
as well the sequence {||/ n || 2 + ||V/„|| 2 }. By the Rcllich-Kondrachov Compactness 
Theorem, see [H Theorem 5.7.1], {/„} has a convergent subsequence in L 2 (T d ). 
From this subsequence, choosing a convergent subsequence of {A„} finishes the 
proof. D 

Lemma 3.4. The image o/I-£ a :J)a^ L 2 (T d ) is dense in L 2 (T d ). 

Proof. By a similar argument to the one found in Lemma 13.11 it is enough to 
show that any smooth function / with support contained in T d \<9A belongs to 
(I — £a)(35a)- Therefore, we need to find a function h in C 2 (T d ) with support in 
T d \<9A such that 

h - Ah = f . 
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From the classical theory of second-order parabolic equations, e.g., see [2j Theorem 
5.7.1], this equation has a smooth solution, which finishes the proof. 

□ 



Proof. (Proof of Theorem |2^|) 

(a) Since 2) a C T-L\, it follows from Lemma |3"TT1 that H\ is dense in L 2 (T d ). 

(b) Denote I— £a = A : T>a —> L 2 (T d ). From Lcmma l3~2l A is linear, symmetric 
and strongly monotone on the Hilbert space L 2 (T d ). By strongly monotone, we 
mean that there exists c > such that 

((AH,H))> c\\H\\ 2 , VHeDa. 

In this case, A satisfy the inequality above with c = 1. By [IT] Theorem 5. 5. a], 
in the conditions above, the Friedrichs extension A : W\ —> L 2 (T 2 ) is self-adjoint, 
bijective and strongly monotone. By an abuse of notation, define now the extension 
C\ : H\ — > L 2 (T 2 ) as (I — A). Since I and A are self-adjoint in W A , this property 
is inherited by C A : U\ -> L 2 (T 2 ). 
For non-positiveness, note that 

((- C A H, H)) = ((-(I - A)H, H)) = -((H, H)) + ((AH, H))>0. 

(c) As mentioned in the proof of (b) above, the Friedrichs extension A : H\ — > 
L 2 (T 2 ) is bijective. So it remains to show that £>a is a core of A : H\ —¥ L 2 (T 2 ). 
For any operator B, denote by G(B) the graphic of B. Then Da is a core for A, if 
the closure of G{A\s a ) l2xl2 in L 2 x L 2 is equal to Q{A). Since A is self-adjoint, A 
is a closed operator, or else, G(A) is a closed set. Thus the closure of G(A\s> A ) is a 
subset of G(A). Let H £ W\, from Lemma [3^41 there exists a sequence {H n } in £>a 



a 



such that AH n converges to AH in L 2 . Hence, as proved in [TTJ Theorem 5.5. 
A^ 1 is a bounded linear operator, and H n converges to H in L 2 , which yields that 
the closure of G{A\® A ) contains G{A). 

(d) Fix a function H in H\ and /i > 0. Put G = (/il — Ca)H. Taking the inner 
product with respect to H on both sides of this equality, we obtain that 

f i((H,H)) + ((-C A H,H)) = ((H,G)) < ((H^H)) 1 ' 2 ((G,G)) 1 ' 2 . 

Since H belongs to 7i\, by (b), the second term on the left hand side is positive. 
Therefore, n\\H\\ < \\G\\ = \\(fil - £\)H\\. 

(e) and (f) We have seen that the operator (I— C\) '■ 2) a —t L 2 (T) is symmetric 
and strongly monotone. By Lemma 13.31 . the embedding TL\ C L 2 (T d ) is compact. 
Therefore, by [TT] Theorem 5.5.c], the Friedrichs extension A : TL\ — > L 2 (T d ), 
satisfies claims (e) and (f) with 1 < Ai < A2 < • • ■ , A„ f °°- I 11 particular, 
the operator —/3a = (A — I) has the same property with < /j,± < /X2 < • • • , 
pi n t 00. Since is an eigenvalue of —£a, a constant function is an eigenfunction 
with eigenvalue 0, then (e) and (f) also hold. □ 

4. Scaling Limit 

Let A4 be the space of positive Radon measures on T d with total mass bounded 
by one endowed with the weak topology. For a measure 7r £ A4 and a measurable 
7r-integrable function H : T d — > R, we denote by (tt, H) the integral of H with 
respect to -k. 
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Recall that {n^ : t > 0} denote a Markov process with state space fijv and 
generator Ln speeded up by TV 2 . Let 7Tj € M. be the empirical measure at time t 
associated to {n^ : t > 0}, which is the random measure in M. given by 

1 



where S u is the Dirac measure concentrated on u. 
Note that 

(7rf.IT> = & J^H(§)^(x), 

for the empirical measures, and (ir, H) = {(p.H)), for absolutely continuous mea- 
sures 7r with L 2 bounded density p, and H E L 2 (T d ). 

Fix T > 0. Let D([0,T],A4) be the space of .M-valued cadlag trajectories 
7r : [0, T] — > M. endowed with the Skorohod topology. Then, the .M-valued process 
{■nf '■ t > 0} is a random element of D([0, T], M) whose distribution is determined 
by the initial distribution of {r]f : t > 0}. For each probability measure p on £In, 
denote by Q^ N the distribution of {vrf : t > 0} on the path space D([0,T],M), 
when t]q has distribution p. 

For a Borel measurable profile 7 : T d — > [0,1], suppose that there exists a 
unique weak solution p of (|2.4p with initial condition 7. We denote by Q^ be 
the probability measure on Z?([0,T],yVf) concentrated on the deterministic path 
ir(t, du) = p(t, u)du. 

Proposition 4.1. Fix a Borel measurable profile 7 : T — > [0,1] and consider a 
sequence {p^ : N > 1} of measures on ttpj associated to 7 in the sense of H2.3\) . 
Then there exists a unique weak solution p of (|2.4p with initial condition 7 and the 
sequence of probability measures Q^ N converges weakly to Q^ as N "f 00. 

It is straightforward to obtain Theorem 12.31 as a corollary of the previous propo- 
sition. The proof of Proposition 14.11 follows directly from the uniqueness of weak 
solutions of (12.41) and the next two results: 



Proposition 4.2. For any sequence {pjq : N > 1} of probability measures with p^ 
concentrated on fijv, the sequence of measures {Q^' N ■ N > 1} is tight. 

Proposition 4.3. Fix a Borel measurable profile 7 : T d — > [0, 1] and consider a 
sequence {pn ■ N > 1} of probability measures on £In associated to 7 in the sense 
of i2.S\) . Then any limit point of Qjjjf is concentrated in absolutely continuous 
trajectories that are weak solutions of (|2.4|) with initial condition 7. 

Proof, (of Proposition \JA\ ). By Proposition 14.21 Q„]f is tight. Thus, if a limit 



point Q* is concentrated in solutions of equation (|2.4[) . by the uniqueness result 
proved in section l4~3l we have that Q* = Q^. And the statement follows from 
Proposition ^. 31 □ 



In Subscction l4.1[ we prove Proposition l4.2l and in Subscction l4.2l we show Propo- 
sition |L3] As a consequence, we have the existence of solutions of (|2.4|) with initial 
condition 7. We complete the proof in subsection (|4.3j) showing the uniqueness of 
weak solutions of (12.41). 
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4.1. Tightness. Here we prove Proposition ^. 21 Let £>([0,T],R) be the space of Un- 
valued cadlag trajectories with domain [0, T] endowed with the Skorohod topology. 
To prove tightness of {k? : < t < T] in D([0,T],M), it is enough to show 
tightness in D([0, T], R) of the real- valued processes {(ir^ , H) : < t < T} for a set 
of functions H : T d — > R which is dense in the space of continuous real functions 
on T d endowed with the uniform topology, see [8]. Futhcrmore, if a sequence of 
distributions in D([0,T],R) endowed with the uniform topology is tight, then it 
is also tight in D([0,T],R) endowed with the Skorohod topology. Here we prove 
tightness of {(n^H) : < t < T} in D([0,T],R), endowed with the uniform 
topology, for H E C 2 {T d ). 

Fix H S C 2 (T d ). By definition {<7rf ,#) : < t < T} is tight in D([0,T],R) 
endowed with the uniform topology if, for the boundedness, 



lim sup P^ 



N 



sup \(tt?,H)\ >m 

0<t<T 



= 0, 



(4.10) 



and, for the equicontinuity, 



Jim limsup P^ N 



sup \(ir?,H)-(ir?,H)\>e 

t-s\<6 



S^Q N 

The limit in (|4.10[) is trivial since 

\(n?,H)\< sup \H(t) 

0<t<T 



= 0, for all e > 0. (4.11) 



So we only need to prove (|4.11[) . By Dynkyn's formula (see appendix in [8]), 

M t N = (Trf , H) - (n», H) - f N 2 L N (^, H)ds (4.12) 

Jo 

is a martingale. By the previous expression, (|4.11[) follows from 



lim limsup 



' flN 



rN 



r Ni 



sup |M t JV -M S JV >e 

\t-s\<5 



, for all e > . 



(4.13) 



and 



lim limsup 



*N 



sup 

0<t-s<8 



N 2 L n (tt{: , H )ds 



<5->0 jv-*oo 

Indeed, wc show the stronger results below: 



> e 



0, for all e > 0. (4.14) 



N 



lim lim sup E, 



sup \M t 

\t-s\<8 



M?\ 



(4.15) 



and 



N 



lim limsup E 



sup 

a<t-s<8 



N z L N {TT j ;,H}ds 



= 0. 



(4.16) 
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To verify (|4.15[) . we use the quadratic variation of MJ* 1 that we denote by (M^). 
By Doob's inequality, we have that 



E 



N 



rN , r N< 



sup \M? - M J a 

\t-s\<5 







r -I 








< 


2E" 

P 1 


sup |M*| 

.0<t<T 










1 


1 


< 


*K» 


sup \M t N \ 2 

_0<t<T 


' < 2<« 


sup (M t N ) 

0<t<T 





Since 

(M t N )= f N 2 [L N {^,H) 2 -2(^,H)L N (n?,H)]ds, 
Jo 

we obtain by a straightforward computation that 
(M t N )= f N 2 J2 J2 ^ + e 1 W^[(Vs(x)-Vs(x + eMH(^)-H(f)) 

Therefore, since ^ x+e . < 1, 



ds . 



T 



< M ?) < w^EE^m^-m) 



i =1 xet%. 



< ™ ( sup \\7H(t) 

N a \ <t<T 



(4.17) 



Thus, Mf converges to zero in L 2 and (|4.15[) holds. 
Wc finish the proof by verifying (|4.16l) . Write 



N 2 L N {^,H) = W ^Y1 E ^ + e ] {{r 1 s{x)-r ls ( X + e ] ))[Hm^)-H{j ! 

= a^E E Vs(x)[e, x+ej (H(^)-H(f))+^ x - e] (H(^)-H(§) 

Define IV C Tfj as the set of vertices whose have some adjacent edge with exchange 
rate not equal to one. Then N 2 Ln(tt^ , H) is equal to 

Iv^E E r, s (x){H(^) + H(^)-2H(jL) 

j=i x(r N 

+lv^E E Vs(x)[^ x+ei (H(^)-H(f))+^ x _ ej (H(^)-H(f)J . 

j=i xer N 

By the Taylor expansion (remember H € C 2 ), the absolute value of the first term 
above is bounded by sup 0<i<T \AH(t)\. Since there are in order of iV^" 1 terms in 
IV, and £ x ,x+et < Ij the absolute value of the second term above is bounded by 

d 

a^E E \H(^)-H(§)\ + \H(^)-H(j-)\<2d sup \\7H(t).e,\. 



j=i xer N 



0<t<T 
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Therefore, there exists C > 0, depending only on H , such that N 2 Ln{it^ , H)\ <C, 
which yields 



N 2 L N (ir£,H)ds 



< C(t - r) . 



and (|4~16|) holds. 



4.2. Caracterization of limit points. Let 7 : T d — > [0, 1] be a Borel measurable 
profile and consider a sequence {pn : N > 1} of measures on f2jv associated to 
7 in the sense of (|2.3|) . We prove Proposition 14.31 in this subsection, i.e., that all 
limit points Q* of the sequence Q^'J^ are concentrated on absolutely continuous 
trajectories ir(t, du) = p(t, u)du, whose density /?(£, u) is a weak solution of the 
hydrodynamic equation (|2.4p with 7 as the initial condition. 



Let Q* be a limit point of the sequence Q^'J^ and assume, without loss of gen- 
erality, that Q^ converges to Q* . 

Since there is at most one particle per site, Q* is concentrated on trajecto- 
ries TTt(du) which are absolutely continuous with respect to the Lebesgue measure, 
TT t (du) = p(t,u)du, and whose density p is non- negative and bounded by 1, se 
Chapter 4]. 

We shall prove the following result: 

Lemma 4.4. Any limit point Q* of Qy^'^ is concentrated is absolutely continuous 
trajectories ir t (du) = p{t,u)du such that, for any H £ 2)a> 



(4.18) 



((pt,H)) -((7, H)) = / ((p s , C K H)) ds . 
Jo 

With the previous lemma we can show Proposition 14.31 



Proof, (of Proposition \4-3\ ). It just remains to extend the equality (|4.18j) to func- 
tions H £ U\. Let H £ H\. Since I - n£ A : £ A ->■ L 2 {T d ) is a core for the 
Fricdrichs extension, there exists a sequence H n £ 3) a such that 

{H n , (I - K £a) #«) -► (#, (I - k £a) H) 

in L 2 (T d ) xi 2 (T d ). Thus, H n ^ H and £ A # n -)■ £ A H, both in L 2 (T d ). Replacing 
H n in equality (|4.18[) . and taking the limit as n — > 00, it finishes the proof. □ 



The remain of this section is devoted to the proof of Lemma 14.41 Fix a function 
H £ 2) a and define the martingale M^ by 

ft 

>ds. (4.19) 



(*?,H) - (*»,H) - f N 2 L N (n?,H), 
Jo 

We claim that, for every 6 > 0, 



lim 1 

N—too 



dA 



sup 

0<t<T 



M, 



N 



> 5 



0. 



(4.20) 



For H £ C 2 , this follows from Chebyshev inequality and the estimates done in the 
proof of tightness, where we have shown that 



lim E 

N— s-oo 



N 



sup \M t 

,0<t<T 



N\ 



< lim E 

AM 



AT 



sup (M t N ) 

0<t<T 



= 0. 



(4.21) 
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For H = h+ A 1a in ®A, the first inequality in ([4.1 7|l is still valid and 



(Mf) < 



A'- 



d 



A 



H(^)-H(f] 



i2 



AT2d-2 



EE [M^)-Mf) 



3=1 a^rV 
d 

' i\/2d-2 / t / J ^x.x+ej -" (. AT j — ■" (at J 



(4.22) 



(4.23) 



where T^v is also defined in the proof of tightness. The expression (|4.22[) goes to 
zero as N increases, since the function h is Lipschitz. For the expression in (|4.23l) . 
let x G T N . If § G A and 2+£i e A c , then ^ +e , < at- The same occurs if f G A c 
and ^jp- G A. If ■&, ^^- both belong to A or A , the exchange rate ^ x+e . is one, 
but \H(^) - F(£)| = \h(^) - h(§)\ < isu Po < t < T |ViJ(i) • e 3 \. In both 
cases, the expression (|4.23|) is of order 0(N~ d ). Therefore, from (|4.21[) . we obtain 

6201). 

The next step is to show that we can replace N 2 1,n by the continuous operator 
C\ in the martingale formula (|4.19|) and that the resulting expression still converges 
to zero in probability. This will follow from the ensuing proposition: 



Proposition 4.5. For any H G 23a, 

ipE \N 2 h N H(f)-C A H(f) 



N-yoo N d 



(4.24) 



xGTi 



Proof. As usual, put H = h+\l\, where h G C 2 (T d ). Rewrite the sum in (|4.24j) 



as 



-L J2 \N 2 L N H(§)-UH(§) +-L J2 \N 2 h N H(j-)-C A H(j-) 



N d 

x<j-T N 

The first term above is equal to 



ier, 



1 E \N>(h(^) + h(^)-2h(fj)-Ah(§] 



N d 



x<£r N 



which converges to zero because heC 2 . The second one is less than or equal to 
the sum of 

w E i a Mat)I 



(4.25) 



xer, 



and 



N d 



nEEk«UW)- fl (*)) 



zerv 3=1 



+mZx- ej (HC-^)-H(§)) 



(4.26) 



Since there are 0(N d x ) terms in Tat, the expression in (|4.25[) converges to zero 
as N — > oo. Since <9A is smooth, the quantity of points x G Tat for which both 



S:r 



A 



x+e. 



and £. 



A 



are different of one is negligible. Therefore, we must only worry 
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about points x G Tn such that, for some j, only one of ^x+e an d £xx-e i s equal 
to iV _1 . This occurs in one of the following four cases: | £ A, ~w L £ -A- and 



x+e 



_ c A c - — g A 



^ G A c and 



;/' + (' 



_ c A • -2L c ,, 



aP, 



*-« g A and H#i G A C : 



JV 



|eA , 2L^1 G A c and ^L^i G A. The analysis of these cases are analogous, thus 
we only consider the first one. Suppose |gA, ^jp~ G A and ^^ G A . In this 
case, the summand in (|4.26[) can be rewritten as 

NC x+e] {H{ x -±^) - #(#)) + N£ x _ e .(H(^) H(§)) 

= l&j ■ e d | [H{ x -±^) - #(§)] + N [H(^) tf (§)] , 
which becomes uniformly (in x G IV) close to 

-\\& d • e,-| sgn(4, • e,-) - ^(f ) - -A&j ■ e, - ^(f ) . 

The condition Vh\g\(u) = — A£(w), which was imposed in the definition of 25a, 
implies that 

lim N£ x+ej {H{^) tf (#)) + JV^ (#(£z£i) - H(f )) = . 

A— »oo 

Therefore, the terms in (|4.26[) converges uniformly to zero, and the same holds for 
the whole sum. □ 



Corollary 4.6. For H G 2) a and for every 5 > 0, 

ft 



lim Qt' N 



sup 

0<t<T 



N 



(*?,H)-(rf,H) 



< , C A H) ds 



o 



> 8 



Proof. By a simple calculation, the martingale defined in (|4.19[) can be rewritten 



a.s 



M t N = (■*?, H) - (ng,H) - [ (ir?,N 2 L N H)ds. 

Jo 
The result follows from proposition 14.51 and expression (|4.20|) . 



D 



From Corollary 14.61 we obtain that (|4.L8p in Lemma T4.4I holds for any function 
H G C{T d ) n ©a- This follows from the fact that for such H, 

ft 



sup 

0<«T 



Jo 



ds 



is a continuous function in D([0,T],A4). 

We still need to show that (|4.18|l holds for any H G Da- For e > 0, define 

(dkf = {u g T d ; dist(u, dA) < e} . 

Let H e be a smooth function which coincides with H in T 2 \(<9A) e and sup T \H e \ < 
sup T | H | . 

Recall that Q* is concentrated on trajectories Wt(du) = p(t,u)du whose density 
p is non- negative and bounded by 1. Then, under Q*, 

sup \{TT t ,H e -H)\ < sup / p(t,u)\H e (u)- H{u)\du 

0<t<T 0<t<TJ(dA) E 

< 2 e d ((dA) £ ) sup \H(u)\. 
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Therefore, for every 8 > 0, 



sup 

0<t<T 



[n t ,H) - (tt ,H) - / (ir s ,£ A H)ds 



> 8 



< 



sup 

0<t<T 



[ir t ,H s ) - (w ,H E ) - / (7T s ,CAH)ds > 8/3 



sup 

0<t<T 



{n u H e -H) > 8/3 



For small enough e, the second probability in the sum above is null. 

If Gi, G2, G3 are continuous functions, the application from D([0,T], M) to 
that associates to a trajectory {7rt,0 < t < T} the number 

rt 



sup 

0<i<T 



(ir t ,Gi) - (7To,G 2 ) - / (tt s ,G 3 } ds 
Jo 



is continuous in the Skorohod metric. Then, 



sup 

0<t<T 



(ir u H s ) - (tt ,H') - / (ir 8) C A H)ds > 8/3 



< lim ( 

W-S-oo 



vA.Af 



sup 

0<t<T 



N TT£\ 



(7Tf,F £ ) - (Tr^iT) 



< , C A H) ds > 8/2 



since Q^'J? converges weakly to Q* and the above set is open. 



By definition, 



sup 

0<t<T 



i^,H e -H) <^J2 \h%x/N)-H{x/N) 



^6Tl 



< 



( £ d ((dA) s ) + O^)) 2 sup \H(u)\, 



because H e coincides with H in T\(<9A) e . Using the same argument as before, we 
obtain 



lim ( 



*A,JV 



sup 

0<i<T 



<^,ff £ > - (7r ,H e ) - f (ir s ,£ A H)ds > 8/3 
Jo 



< lim 

A-s-oc 



nA.iV 



sup 

0<t<T 



(n t ,H) - (n ,H) 



(tt s ,CaH) ds > 8/9 



+2 lim 



l,A,A 



N- 



sup 

0<t<T 



[<K t ,H s -H) > 8/9 



Again, for small enough e, the second probability in the sum above is null. Re- 
calling the corollary 14.61 we finally conclude that Q* is concentrated on absolutely 
continuous paths nt(du) = p(t,u)du with positive density bounded by 1, and Q* 
a.s. 

((pt,H)) - ((po, H)) = f {{p. , C A H)) ds , 
Jo 

for any H £ 2) a. Therefore we have proved Lemma WM 
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4.3. Uniqueness of weak solutions. Now, we prove that the solution of (|2.4|) is 
unique. It suffices to check that the only solution of (|2.4[) with po = is p = 0, 
because of the linearity of C\. Let p : K + x T d — > R be a weak solution of the 
parabolic differential equation 

d t p = Cap 
P(0,-) = 0. 



By definition, 

({p t ,H)) = I ((p s ,C A H))ds, (4.27) 



'o 

for all functions H in T~L\ and all t > 0. From the theorem 12. 2[ the operator — C\ 
has countable eigenvalues {/i„ : n > 0} and eigenvectors {F n }. All eigenvalues 
have finite multiplicity, = /xq < pi < ' • ' j an d hm„_ ! , 00 p n = oo. Besides, the 
eigenvectors {F n } form a complete orthonormal system in the L 2 (T d ). Define 

neN n u + /^ 

for all £ > 0. Notice that -R(O) = and R(t) is well defined because p t belongs to 
L 2 {T d ). Since p satisfy (J437) . we have that f t ((p t , F„)) 2 = -2p n ((p t ,F n )) 2 . Then 

because J2 n < N ^T+jt) (&*> F ^)) 2 converges uniformly to J^neti n^i+^ dPt^^)) 2 ^ 
when N increases to infinity. Thus R(t) > and {-^R)(t) < 0, for all t > and 
R(0) = 0. From this, we obtain R(t) = for all t > 0. Since {i 7 ^} is a complete 
orthonormal system, (pt, pt} = 0, for all t > 0, which implies p = 0. 

Acknowledgments: The authors T. Franco and A. Neumann would like to thank 
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